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ABSTRACT 

We present measurements of two types of cluster galaxy alignments based on a volume limited and 
highly pure (> 90%) sample of clusters from the GMBCG catalog derived from SDSS DR7. We detect 
a clear BCG alignment (the alignment of major axis of the BCG toward the distribution of cluster 
satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG 
absolute magnitude decrease, and becomes weaker as BCG stellar mass decreases. No dependence 
of the BCG alignment on cluster richness is found. We can detect a statistically significant (> 3 
sigma) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the 
BCG) only when we use the isophotal fit position angles (PAs, hereafter), and the satellite alignment 
depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests 
the detected satellite alignment based on isophotoal PAs from the SDSS pipeline is possibly due to 
the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply 
interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with 
our current photometric SDSS data. We perform our measurements on both SDSS r band and i band 
data, but did not observe a passband dependence of the alignments. 

Subject headings: Galaxies: clusters: general - large-scale structure of universe 



1. INTRODUCTION 



Galaxy orientations contain important information 
about the gra vitational environment in which it resides. 
Brown (19381) pointed out that galaxy orientations may 
not be iso tropic due to the large s cale gravitational in- 



teraction. |Hawley fc Peebles ( 1975 ) reported a weak ev- 
idence of anisotropy of the galaxy orientation. Galaxy 
orientation becomes especially interesting in the vicinity 
of galaxy clusters, where the strong gravitational field 
may produce detectable orientation preference for both 
central galaxies and satellite galaxies. With the advent 
of mode rn sky surveys, su ch as Sloan Digital Sky Survey 
(SDSS) ( |York et al.||2000| , the shape and orientation of 
galaxies can be measured to high precision. The large sky 
coverage substantially increases the sample size to allow 
statistically significant measurements on cluster galaxy 
alignments, from which our understanding of the cluster 
formation process can be greatly improved. 

The term galaxy alignment has been used extensively 
in the literature and refers to alignments in different con- 
texts. In the galaxy cluster environment, there are two 
types of alignments that are of great interest. The first 
is the alignment of the major axes of the cluster satel- 
lite galaxies towards the cluster center, which we will call 
"satellite alignment" . In our case, as an operational defi- 
nition, we will consider the alignment between the major 
axes of the satellite galaxies and the BCG 5 . The second 
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type of alignment is the alignment of the BCGs major 
axis towards the distribution of satellite galaxies in the 
cluster. We will call this alignment "BCG alignment" 
hereafter. 

In addition to these two types of alignments, the possi- 
ble alignments between t he major axises of th e satallitc 
galaxies and the cluster (Plionis et al. 2003[), and be- 



tween the cluster shape and large scale structures ( Paz 
eF^[2008l IFaltenbacher et al |[2009l |Wang et al||20tM 
Paz et al. 120111) have been studied, but these are beyond 
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the scope ot this current paper. 

There are extensive studies on these two types of align- 
ments with both simulations and observations. It is ar- 
gued, based on simulations, that the preferred accretion 
direction of satellite halos toward the host halo along 
the fil aments is largely responsible for the BCG align- 
ment (|Tormen || 1997 I Vitvitska et al. || 2002l IKnebe et al 



2004| |Zentner eTaTf2005| |Wang et al.||2005 

tions of this align ment from rea l data has been reported 
by many teams (Sastr 

T 



Dressier 
Brainer 



Carter 



2 



19G8 

Metca 



2006 



1978] 

2005l|Vangetal. 
et al.||2008| |Siverd et al.||2009 

2010), thouglTnon-detection o: 



he detec- 
P 1 

Austin & Peach} [1974 
i'el \iM\ iB inggeli 1982 



lAzzaroet al.|2 007 Wang 
|Niederste-Ostholt et al. 



reported ( |Tucker & Peterson| 1988| |Ulmer et al.||1989[ ) 



this alignment were also 



play a major role in its formation ( Ciotti & Dutta 
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DSS data has been re- 



especially optical data, is not unambiguous. The central galaxy in 
a cluster (the one which resides near the bottom of the cluster po- 
tential well) is very often the brightest galaxy (BCG) in the cluster. 
This BCG is then coincident with the region with the deepest po- 
tential traditionally identified in theory as the center of a cluster. 
Using the BCG as the cluster center simplifies precise comparisons 
between observations and theory. 
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2009 ). In Table [Tj we summarize the previous work 



that reports the existence and non-existence of these two 
types of alignments based on real data. On the other 
hand, the intrinsic alignment of galaxies will contaminate 
gravitational lensing measurements and therefore needs 
to be carefully modeled in lensi n g analysis. A l ong th ese 
lines, |Mandelbaum et al.| ( |2006[ ) ; |Hirata et aL] ( |2007[ ) re- 
ported correlations between intrinsic shear and t he den- 



sity field ba sed on data from SDSS and 2SLAQ (|Croom 
eTaLl[2009 ). 



in general, the cluster galaxy alignment signals are 
weak, and their measurement requires high quality pho- 
tometry and well measured galaxy PAs. Moreover, a 
galaxy cluster catalog with high purity, well-determined 
BCGs and satellite galaxies is important too. Since mea- 
suring the redshift evolution of the alignment is crucial 
for understanding its origin, the cluster catalog needs 
to be volume limited and maintain constant purity for 
a wide redshift range. Most of the existing alignment 
measurements (see Table [lj are based on galaxy clus- 
ters/groups selected from spectroscopic data. In SDSS 
data, due to the high cost of obtaining spectra for a 
large population of galaxies, the completeness of spec- 
troscopic coverage is limited to r band Petrosian mag- 
nitude r r <17.7, corresp onding to a median redshift of 



0.1 (Strauss et al. 2002). This greatly limits the ability 
to look at the redshift evolution of the alignments. 

On the other hand, one can also measure the align- 
ments by using photometrically selected clusters. The 
advantage of photometrically selected clusters lies in the 
large data sample as well as relatively deep redshift cov- 
erage, allowing a study on the redshift evolution of the 
alignments. However, there are clear disadvantages too. 
For example, the satellite galaxies are prone to contami- 
nation from the projected field galaxies, which will dilute 
the alignment signals. The level of this contamination 
may also vary as redshift changes, complicating the in- 
terpretation of the alignment evolution. 

In this paper, we show our measurements of the two 
types of alignments based on a volume limited and highly 
pure (> 90%) subsample of clu sters from the GM BCG 
cluster catalog for SDSS DR7 ( |Hao et al.||2010[ ). The 
large sample of clusters allows us to examine the depen- 
dence of the alignments on cluster richness and redshift 
with sufficient statistics. With this catalog, we detect a 
BCG alignment that depends on redshift and the abso- 
lute magnitude of BCG, but not on the cluster richness. 
We also observe that the satellite alignment depends on 
the apparent brightness of the BCG and the methods 
the position angles are measured. We can only see a 
statistically significant satellite alignment at low redshift 
when we use the isophotal fit PAs (see § 3.3 for more 
details). Furthermore, we notice that the satellite align- 
ment based on isophotal fit PAs depends strongly on the 
apparent magnitude rather than the absolute magnitude 
of the BCG. This suggests that the measured satellite 
alignment is more likely due to the isophotal fit PAs, 
whose measurements are prone to contaminations from 
the diffuse light of the BCG. 



The paper is organized as following: in §2, we introduce 
the two parameters used to quantify the two types of 
alignments. In § 3, we introduce the data used in this 
paper. In § 4, we present and discuss our measurement 
results. By convention, we use a ACDM cosmology with 
h = 1.0, £1,„ = 0.3 and Q\ = 0.7 throughout this paper. 
All angles measurements that appear in this paper are 
in units of degrees. 

2. ALIGNMENT PARAMETERS 

In this paper, we consider two types of alignments: (1) 
Satellite alignment; (2) BCG alignment. Each of them 
is quantified by a corresponding align ment parameter. 
For the sa t ellite alignment, we f ollow |Struble &: Pee-| 
pes|(IT985l);|Pereira & Kuhn (p005l) and use the following 
alignment parameter: 



5 = 



N 



45 



(1) 



where <pi is the angle between the major axes of the satel- 
lite galaxies and the lines connecting their centers to the 
BCGs, as illustrated in the left panel of Figure [TJ N 
is the number of satellite galaxies in the cluster. For 
every cluster, there will be a unique alignment parame- 
ter S measured, which is the mean angle (j) of all cluster 
satellite galaxies subtracted by 45. If the major axes of 
satellite galaxies do not preferentially point to the BCG 
of the cluster, the 4> will randomly distribute between 
-45 and 45 (degrees), leading to S = 0. On the other 
hand, if the major axes of satellite galaxies preferentially 
point to the BCG, there will be 5 < 0. The standard 
deviation (error bar) of 5 can be naturally calculated by 

Serr = ^=1 {&~5~Ahf/N. 

For the BCG alignment, we will focus on the angle 
9 between the BCGs PA and the lines connecting the 
BCG to each satellite galaxy, as illustrated in the right 
panel of Figure [JJ In our cluster sample, each cluster has 
more than 15 satellite galaxies (see §3.2 for more details). 
Therefore, instead of looking at the full distribution of 9 
from all clusters, we will focus on the mean of 9 measured 
for each cluster. In analogy to the satellite alignment 
parameter 5, we introduce a BCG alignment parameter 
7 defined as: 
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where 9i is the angle between BCGs PA and the line 
connecting the BCG to the i fh satellite galaxy (see the 
right panel of Figure [lj). N is the number of satellite 
galaxies in the cluster, similarly, each cluster will corre- 
spond to a BCG alignment 7, which is the mean of angle 
9 subtracted by 45. If the major axis of BCG preferen- 
tially aligns with the distribution of the majority of the 
satellite galaxies, the 7 < 0. If no such preference exist, 
7 = will be expected. The uncertainty of 7 can be 

readily calculated as j err = i®i — 7 ~~ 45) 2 /7V. 

These two parameters quantify the two types of align- 
ments and are easy to measure. In the follows, we will 
focus on these two quantities and their dependencies on 
various cluster/BCG properties. 
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Figure 1. left: Illustration of the angle <j> used in the definition of satellite alignment parameter <5; right: Illustration of the angle 9 used 
in the definition of BCG alignment parameter 7. 



3. DATA 

In this section, we describe the details of the galaxies, 
their PA measurements, and the galaxy cluster sample 
used in our measurement. 



3.1. Galaxy Catalog 

The galaxies we use are from the Data Release 7 of the 
Sloan Digital Sky Survey ([York et al.|2000l JAbazaj ian & 
Sloan Digital Sky Survey||2008p . The SPSS is a multi- 
color digital CCD imaging and spectroscopic sky survey, 
utilizing a dedicated 2. 5- meter telescope at Apache Point 
Observatory, New Mexico. It has recently completed 
mapping over one quarter of the sky in u,g,r,i and z 
filters. DR7 is a mark of the completion of the original 
goals of the SDSS and the end of the phase known as 
SDSS-II. It includes a total imaging area of 11663 square 
degrees with 357 million unique objects identified. 

This work focuses on the Legacy Survey area of 
SDSS DR7, which covers more than 7,500 square de- 
grees of the North Galactic Cap, and three stripes 
in the South Galactic Cap totaling 740 square de- 



grees (Abazajian & Sloan Digital Sky Survey 2008). 
The galaxies are selected from the PhotoPrimary view 
of the SDSS Catalog Archive Server with object type 
tag set to 3 (galaxy) and i-band magnitude less than 
21.0. Moreover, we require that the galaxies did 
not trigger the following error flags: SATURATED, 
SATUR_CENTER, BRIGHT, AMOMENTJVIAXITER, 
AMOMENTJ3HIFT and AMOMENTTAINT. 

In addition to the above selection criteria, we also re- 
ject those galaxies with photometric errors in r and i 
band greater than 10 percent. Additionally, we require 
the ellipticity 6 of each galaxy in the r-band and i-band 
to be less than 0.8 in order to remove edge-on galaxies 
whose colors are not well measured. By doing this, we 
will retain about 95% of the total galaxies. All the mag- 
nitudes used in thi s paper are dust extinction corrected 
model magnitudes (Abazajian & Sloan Digital Sky Sur- 
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3.2. Galaxy Clusters 

In order to measure the two types of alignments, we 
need to have a galaxy cluster catalog with well deter- 
mined member galaxies. The GMBCG cluster catalog 
for SDSS DR7 is a large catalog of optically selected 
clusters from SDSS DR7 using the GMBCG algorithm. 
The catalog is constructed by detecting the BCG plus 
red sequence feature that exists among most clusters. 
The cluster satellite galaxies are within 2a of the red 
sequence mean color detected using a Gaussian Mixture 
Model. Since the red sequence of each cluster is measured 
individually, it allows a more accurate satellite galaxy se- 
lection than using a universal red sequence model. We 
count the satellite galaxies down to 0.4L* at the cluster's 
redshift. In Figure [2] we show the color distribution of 
galaxies around a cluster in GMBCG catalog. To study 
the redshift evolution of the alignment, we choose the vol- 
ume limited sample of clusters with redshift below 0.4, 
where the satellite galaxies are selected using g — r color. 
Also, to get high purity, we choose clusters with rich- 
ness equal to or greater than 15, which leads to a purity 
above 90% across the redshift range. As a result, there 
are about 11,000 clusters with over 260,000 associated 
satellite gal axies. More detai ls of the cluster catalog can 
be found in Hao et al. (2010). 



3.3. PA Measurements 

In the SDSS data reduction pipeline, the PAs of 
gala xies are measured wi th several different meth- 
ods (Stoughton et al. 2002). In this work, we will use 
three of them: isophotal PA, exponential fit PA 7 and 
De Vaucouleurs fit PA. 

To measure the isophotal PA, the SDSS pipeline mea- 
sures out to the 25 magnitudes per square arc-second 
isophote (in all bands). The radius of a particular 
isophote as a function of angle is measured and Fourier 
expanded. From the coefficients, PA (isoPhi) together 
with the centroid (isoRowC,isoColC), the major and mi- 
nor axes (isoA,isoB) are extracted. 

For the exponential fit PA, the SDSS pipeline fit the 
intensity of galaxy by 



I(R) =Ioexp[-lM(R/Reff) 



(3) 



where the profile is truncated outside of 3R e ff and 
smoothly decreases to zero at 4R e ff. After correcting 
for the PSF, the PA is calculated from this fitting and 
reported as expPhi in the CASJOB database. The De 
Vaucouleurs fit PA follows the same procedure as expo- 
nential fit PA, except the model in the fitting is 



I(R) =I exp[-7.67{R/R eff ) 



l/4i 



(4) 



where the profile is truncated outside of 7R e ft and 
smoothly decreases to zero at 8i? e //. The PA from this 
fitting is reported as devPhi. For more details about 



these P A measurements, one can refer to Stoughton et al. 
(2002). All the angles are in degrees East of North by 



7 Since most satellite galaixes are selected using red sequence, 
the exponential fit may not be suitable. In this paper, we want 
to demonstrate how will the PA measurement affect the alignment 
signal, and therefore include it in our discussion. 



convention. By comparing these different PAs, isophotal 
PA tends to trace the exterior shape of the galaxy while 
the two model fit PAs tend to trace the inner profile of 
the galaxy. At low redshift, the BCG is very bright and 
its diffuse light may severely affect the measurement of 
the outer part of the nearby galaxies. Therefore, the 
isophotal PA is more susceptible to this artifact, leading 
to an "artificial" orientation preference toward BCG. To 
give a sense of the BCG diffuse light effect, we show a 
low redshift and a higher redshift cluster in Figure [3] 

3.4. Control Sample 

As we talk about the detection of the alignment sig- 
nals, we need to have a control sample to compare with. 
Naively, one may directly compare the measured align- 
ment signal to what is expected for non-detection. How- 
ever, this will implicitly assume that the data as well as 
the measurements are free from any systematics. When 
talking about the significance level of the detection, it is 
not sufficient to consider only the statistical uncertain- 
ties. Systematics (intrinsic scatter), if they exist, also 
need to be considered. Therefore, introducing an appro- 
priate control sample and applying the same measure- 
ments we apply to the cluster data will help to eliminate 
the possible false detection resulted from potential sys- 
tematics. 

For this purpose, we prepare our control sample as fol- 
lows: we shuffle the BCGs by assigning random positions 
(RA and DEC) to them, but keep all other information 
of the BCGs unchanged. Then around each BCG (at 
new random position), we re-assign "cluster satellites" by 
choosing those galaxies that are falling within the R sca ie 
from the BCG. Also, the i band magnitude of the "satel- 
lites" should be in the range from 14 to 20. The R SC aie 



measured i n Mpc, plays the role of virial radius (Hao 



et al. 2010). In addition to this non-cluster sample, we 
also add another random PA to every galaxy in both the 
true cluster sample and the non-cluster sample by re- 
placing the measured galaxy's PA with a random angle 
uniformly sampled between and 180 degree. This ran- 
dom PA will serve as another control sample to double 
check the possible systematics in our measurements. 

4. RESULTS 

4.1. Satellite Alignment 

There are two basic questions concerning the forma- 
tion of the satellite alignment: (1) is it a residual feature 
of the initial condition of the cluster formation? (2) or is 
it a dyn amically evolving effect that varies as the cluster 
evolves ( |Pereira fc K uhn 2005), for example, due to the 
tidal torque? The two different scenarios lead to differ- 
ent redshift dependence of S. If it is left over from the 
initial alignment, its strength should decrease as redshift 
decreases. On the other hand, we should see a stronger 
alignment sign al at low redshift if it is a dynamically 



evolving effect (Ciotti & Dutta 1994 Catelan et al. 2001 
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Kuhlen et al. 2007; 
2010). Therefore, L 



Pereira et al.||2008[ |Pereira Bryan| 
ooking at the redshift dependence of 
the measured 6 is our primary interest. To do this, we 
first measure the S for each cluster using 4 different PAs, 
i.e., the random PA, exponential fit PA, De Vaucouleurs 
fit PA and isophotal PA in r band. We bin the clusters 
into redshift bins of size 0.05. In each bin, we calculate 



Intrinsic alignment 



5 



ObjlD: 58B0171 
yu: 1.884 1. 
a: 0.102 0.516 
photoz: 0.259 




slope:-0.060 
intercept: 3.009 



17 



■ a 

mag 



' 9 



?0 



Figure 2. ie/i; Galaxy g— r color distribution around a cluster overlaid with a model constructed of a mixture of two Gaussian distributions. 
The red curve corresponds to the red sequence component while the blue one corresponds to the sum of background galaxies and blue 
cluster satellites. The green vertical line indicates the color of the BCG. /i and <r are the means and standard deviations of the two Gaussian 
components, right: Color-magnitude relation for the same galaxies. Galaxies within the 2cr clip of the red sequence component are shown 
with red points; the green line indicates the best fit slope and intercept of this red seqence. 
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Figure 3. Left is a cluster at redshift 0.103 and the right is a cluster at redshift 0.310. The light from BCG will pose a risk on the proper 
measurements of the PAs of the satellite galaxies, especially at low redshift. 



the weighted mean of 5 and the standard deviation of 
the weighted mean, with the weights specified by 1/ S^ rr . 
Then, we perform the same measurement on the control 
sample, and present the results in Figure [4] and Figure [5] 
Based on the results in Figure |4j one can see that the d 
measured using isophotal PAs deviated from that based 
on random PAs in a statistically significant way at low 
redshift, though approaching zero as redshift increases. 
While the 5 measured using the exponential fit PA and 
De Vaucouleurs fit PA are consistent with that measured 
using random PAs exce pt in the lowest redshift bin. In 
Pereira & Kuhn (2005), the authors used isophotal fit 
PAs and also considered a cluster sample (sample B) 
with satellite galaxies selected using red sequence. When 
limiting the r band magnitude to less than 18, they mea- 
sured a S = —1.06 ± 0.37, which is consistent with our 
results in the lowest redshift bin using isophotal fit PAs. 
The results based on the random control sample in Fig- 
ure [5] show that S measured using all types of PAs are 
consistently zero across the redshift range. In Figure [6j 
we plot the axis ratio of the satellite galaxies in different 
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Figure 4. Satellite alignment for clusters at different redshift bin 
of size 0.05. The legend random PA indicates using of randomized 
PAs. exp indicates using exponential fit PAs. dev indicates using 
De Vaucouleurs fit PAs and iso indicates using isophotal PAs. r in- 
dicates the SDSS r filter. This legend convention is also applicable 
to other figures in the paper. 
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trol sample. 
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Figure 7. The dependence of satellite alignment 8 on the r-band 
apparent magnitudes of the corresponding BCGs. This shows a 
strong dependence of 8 measured using isophotal PAs on the ap- 
parent magnitude of the BCGs. 
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Figure 6. The axis ratio of the satellite galaxies. Here, the error 
bar is the standard deviation to the mean in that bin. The legend 
satellite expAB_r refers to the b/a ratio is measured by fitting the 
satellite galaxy with an exponential profile in r band. devAB_r 
and isoAB_r refer to the b/a ratio by fitting De Vaucouleurs and 
isophotal profiles respectively. This convention is also applicable 
to other figures in this paper. 



redshift bins. The axis ratio determines how precisely the 
PAs are measured. From Figure [6j the axis ratio mea- 
sured using the isophotal method aoes not vary much as 
redshift increases, indicating that the diminishing satel- 
lite alignment based on isophotal fit PAs is not due to 
the decreasing S/N at high redshift. 

There are two possible explanations for the measured 
6 in the low redshift bins. The first one is that the dif- 
fuse light from the BCGs affects the measurements of 
the PAs of the cluster satellite galaxies'. This creates an 
artificial preference of major axes of the satellite galax- 
ies. This contamination is most severe when the PAs are 
measured using isophotal fit, but less prominent when 
the PAs are measured using exponential fit and De Vau- 
couleurs fit. This is because the isophotal PAs are sensi- 
tive to the shape of the outer profile of galaxy while the 
model fit PAs are more determined by the inner profile 
of the galaxy. 

The second possible explanation to these results is the 
twisting of galaxy. This leads to different PAs when we 
use different methods. The outer rim of the galaxy is 
more susceptible to the tidal torque so that the align- 
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Figure 8. The dependence of satellite alignment 8 on the r-band 
absolute magnitude of the BCGs. This shows that there is very 
little dependence of 8 measured using isophotal PAs on the r band 
absolute magnitude of the BCGs 

ment will show up when we use the isophotal fit PAs. 
One way to distinguish these two explanations is to look 
at the way S depends on the absolute and apparent mag- 
nitudes of the BCG. To see this, we plot the measured S 
for all the clusters with respect to apparent magnitudes 
and absolute magnitudes of the corresponding BCGs in 
the Figure [7] and Figure [8] respectively. One can see that 
S shows a strong dependence on the apparent magnitude 
but not on the absolute magnitude. Therefore, we con- 
clude that the 8 in the low redshift bins in Figure [4] is 
more likely resulted from the artifact of the PA measure- 
ment. In Figure [9l we show the absolute magnitude vs 
redsfhit for the BCJGs. 

4.2. BCG Alignment 
4.2.1. Redshift Dependence 

Now, we consider the BCG alignment. We first look 
at the redshift dependence of 7. Again, we perform our 
measurements on both the cluster sample and the ran- 
dom control sample. The results are shown in Figure [TO] 
and Figure [TT1 

The results from the random control sample show that 
7 is consistently zero. From the cluster sample, we de- 
tect a clear BCG alignment and its strength decreases as 
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Figure 9. Photometric redshift vs. r-band absolute magnitude 
for BCGs. The red overplotted dots are the means in each redshift 
bin of size 0.05. 
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Figure 12. The axis ratio b/a of BCGs and clusters vs. redshift. 
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Figure 10. BCG alignment at redshift bins of size 0.05 using 
different PAs measured from the cluster sample. 
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Figure 11. BCG alignment at redshift bins of size 0.05 using 
different PAs measured from the random sample. 



redshift increases. From the figures, we can also see that 
7 is almost the same no matter how the PAs of BCGs 
are measured. This has two important implications: 1. 
the PAs of BCGs are well measured by both isophotal fit 
and model fit; 2. the diffuse light of satellite galaxies does 
not affect the PA measurement of the BCG. Before we 
can conclude the redshift dependence of 7, we still need 
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Cluster PA 
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Figure 13. The distribution of BCGs' and clusters' PAs. 

to do one more test. That is, if the axis ratio (b/a) of 
the BCG and cluster become systematically smaller due 
to the decreased S/N at higher redshift, the measured 
strength of 7 will decrease too. 
To calcula te the clus t er b/ a, we use the method as 



describ ed in |Kim et al.| ([2002]); |Niederste-Ostholt et al. 
(2010). For clusters, the axis ratio is defined as b/a = 

v/Q 2 + U 2 )/(l + v/Q 2 + U 2 ), with the stokes pa- 



U 2 ) 
and U 



2M T 



The radius- 

= (x 2 /r 2 ), 



(1 

rameters Q = M xx — M Va 
weighted second moments are given by M a 
M yy = (y 2 /r 2 ) and M xy — (xy/r 2 ) with r 2 = x 2 + y 2 . 
x and y are the distances between the satellite galaxies 
and BCG in the tangent plane. For BCGs, we use the 
measured b/a in SDSS pipeline based on isophotal fit, 
exponential fit and DeVoucular fit. In Figure [l2j we plot 
the b/a in each redshift bin of size 0.05. From the results, 
the axis ratio does not depend on redshift in a statisti- 
cally significant way. Furthermore, we also checked that 
both the BCGs' PAs and cluster PAs are distributed ran- 
domly, as show in the right panel of Figure[l3] Therefore, 
the evolution of 7 in Figure [10] should not result from the 
S/N variation of the PA measurements. 

4.2.2. Magnitude Dependence 

We measure the BCG alignment vs the r band absolute 
magnitudes of the BCG. The results are presented in 
Figure [T4| 

Fronithe plot, we see that 7 strongly depends on the 
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Figure 14. The dependence of BCG alignment on the absolute r 
band magnitudes of BCGs. 
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Figure 15. The axis ratio b/a of BCGs and clusters vs. the r 
band absolute magnitude of BCG. 



BCGs absolute magnitude. To further show that this is 
not due to S /N of the measurement of BCGs shape, we 
plot the axis ratio b/a of BCGs and clusters as a func tion 
of the BCGs r band absolute magnitude in Figure p~5| 
There is no dependence of b/a on the BCGs absolute 
magnitude. 

From the above results, we see that the 7 depends on 
both redshift and the BCG absolute magnitude. To show 
this more clearly, we bin the cluster samples into photoz 
bins of size 0.05 and absolute magnitude bins of size 0.5. 
Then, we calculate the mean 7 in each bin and plot the 
results in Figure [T6| The color in the plot indicates 7. 

From the plot, we can see that the 7 increases (i.e. ab- 
solute signal decreases) as redshift increases and as the 
absolute magnitude of the BCG decreases. Since the ab- 
solute magnitude of a galaxy is proportionally correlated 
with its mass, the above results indicate that the more 
massive BCGs tend to be more aligned with the cluster 
orientation. A subsample of the BCGs (~2800 BCGs) 
has their stellar masses measured in the MPA-JHU value- 



2003 



added catalog fo r SDSS DR7 ( Kauffmann et al., 
Salim et al.||2007 ). This allows us to directly look at the 
trend ot 7 with respect to the BCG stellar mass (a proxy 
of the total mass) and redshift. We choose three stellar 
mass bins and plot the 7 vs. redshift in each of them in 
Figure 17 We can see that the 7 increases as redshift 



shows more negative 7. 

4.2.3. Richness Dependence 

Next, we check whether 7 depends on cluster richness. 
To do this, we bin the clusters by their richness in bins 
with edges at 15, 25, 35, 50, 65 and 90. We did not go 
to the richness bins above 90 due to the smaller number 
of clusters in that richness range. Then, we look at the 
mean 7 in each bin. Note that we did not re-bin them 
into different redshift bins to keep the number of clusters 
reasonably large. This will not affect our purpose for de- 
tecting richness dependence since clusters from different 
redshifts are randomly falling into different richness bins. 
We plot the results in Figure l8l Based on the results, 



increases in each stellar mass bin and the massive BCGs 



we do not see a dependence on cluster richness. 

4.3. Redshift Evolution of BCG Alignment Once Again 

Measuring the redshift evolution of BCG alignment is 
very important for understanding its origin. However, 
the redshift evolution of the measured alignment signal 
needs to be interpreted with great caution, especially 
for the cluster samples selected using photometric data. 
There are at least four factors that will introduce system- 
atic redshift dependence and complicate the interpreta- 
tion. (1) The S/N of the galaxy shape measurements 
will decrease as redshift increases. (2) When we look at 
clusters of different redshift, we need to make sure we 
are comparing the same population of satellite galaxies. 
That is, the cluster catalog needs to be volume limited in 
the redshift range. (3) The purity of clusters need to be 
consistently high across the redshift range. The change 
of purity will lead to decreased mean alignment signal. 
(4) The level of contamination from the projected field 
galaxies are prone to redshift dependence. Higher level 
of contamination will dilute the alignment signal. We 
have addressed (1) and (2) in previous sections, where 
we introduce the results. In the follows, we will focus on 
the (3) and (4). 

The alignment signal from the falsely detected clus- 
ters should be consistent with zero, decreasing the mean 
alignment of the whole sample. For the subsample of 
GMBCG clusters with richness equal or greater than 15, 
it has been shown that the purity is consistently above 
90% and does not vary more than 10% across the redshift 
range from 0.1 to 0.4 ( |Hao et al.|2010| . To further show 
the purity variation will not produce the observed red- 
shift dependence of 7, we choose another subsample of 
clusters with even higher purity. We choose clusters with 
richness greater than 25, which have a purity of above 
95% and vary less than 5% in the redshift range. In Fig- 
ure |19[ we plot the BCG alignment parameter 7 from 
this subsample. Though the overall signal level increases 
a little at low redshift, the trend of redshift dependence 
does not differ much from the full sample with a lower 
richness threshold 15. So, the purity change should not 
explain the measured strong redshift dependence of 7. 

On the other hand, satellite galaxies selected using 
red sequence colors have different levels of contamination 
from projected field galaxies as redshift changes. This is 
mainly caused by the different degree of overlap between 
the red sequence population and the fie l d gala xy popu- 
lation (e.g. see Figure 14 in Hao et al. (2010)). There 
are two competing effects that will increase or decrease 
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Figure 17. Redshift dependence of 7 in three BCG stellar mass bins. 
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Figure 18. Dependence of BCG alignment on cluster richness. 
Here, we bin the richness into bins with edges at 15, 25, 35, 50, 65 
and 90. The results do not show a statistically significant depen- 
dence of BCG alignment on cluster richness. 

the contamination. First, the separation between the red 
sequence component and the field galaxy component. As 
redshift increases, the two components separate farther, 
leading to decreased projection contamination in the red 
sequence. The second effect is the broadening of the dis- 
tribution of both red sequence and field galaxy. As red- 
shift increases, the measured width of red se quence in- 
creas es mainly due to the photometric errors ( |Hao et al.| 
|2009| . This will increase the chance of projected field 
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Figure 19. BCG alignment at redshift bins of size 0.05 based on a 
subsample with higher purity and lower purity variations across the 
redshift range. The results indicate that the redshift dependence 
are not affected by the purity of the cluster sample we are using. 

galaxies being identified as satellites when we select the 
satellite galaxies by color. Therefore, the actual contam- 
ination level is the compromise of these two effects. 

We can describe the measured alignment parameters, 
7 8 , as a combination of alignment from real cluster satel- 

8 Since we did not see a significant <5 signal, we will consider only 
7 in all the discussions hereafter. But the method described can 
also be applied to the 8 case. 
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lites and projected field galaxies. If we denote the align- 
ment parameters from our measurements as j m , then we 
can decompose it into two parts as follows: 



7rc 



N c + N f 



45 



(5) 



where N c is the number of true cluster satellite galax- 
ies and Nf is the number of projected field galaxies. 
We can introduce the fraction of real cluster satellite 
as f c ( z ) = N C /(N C + N f ), the BCG alignment from 

true cluster members as j c = l^c — 45 and 

the BCG alignment from the projected held galaxies as 
7/ = J2^=a Qj/Nf — 45. Substitute these definitions into 
Equation [5] and take ensemble average of the clusters, we 
will have: 



(7™) = (fc(z)) ( 7 c> + [1 - (f c (z))] (7/) 



(6) 



where (...) denotes the average over the cluster ensemble. 
As the mean alignment signal from the field is consis- 
tent with zero, the alignment parameter 7 from the true 
cluster satellites is related to the measured one through 
the redshift dependent fraction f c (z). To the first order 
approximation, we can separate f c (z) into two parts as 
fc{z) = fconst x f(z), where f con st is a redshift indepen- 
dent component of the fraction, indicating the "intrinsic" 
fraction of true satellite based on color selection. f(z) is 
the redshift dependent part, corresponding to the effect 
we described above. Then, the redshift dependence of 
the measured alignment 7 will be mainly determined by 

/(*)• 

In the GMBCG catalog, we also measured a weighted 
richness, which takes into account the different degree 
of overlaps between red sequence and th e field galax- 



ies at different redshift (Hao et al. 
ence between weighted richness anc 



2010]). The diffcr- 
the direct member 



count richness is a good estimator of the number of pro- 
jected galaxies due to the effect described above. The 
fraction of contamination can therefore be estimated by 
the ratio of this di ffere nce to the direct member count 
richness. In Figure [20j we plot the fraction of contam- 
ination (1 — (f(z))) as a function of redshift in bins of 
size 0.05. The fraction is almost constant except for the 
lowest redshift bin. Again, this cannot explain away the 
dependence of 7 on redshift as shown in Figure [10] and 
Figure [l9j Therefore, after considering all the possible 
systematics known to us, the measur ed redshift depen- 
dence of 7 st ill cannot be explained. In Nicdcrstc-Ostholt 
et al. ( 2010[ ), the authors also reported a different BCG 
alignment between one low redshift bin (0.08 - 0.26) and 
another high redshift bin (0.26 - 0.44), which is consistent 
with the results we find here. 

4.4. Conclusions and Discussions 

We measure the satellite alignment and BCG align- 
ment based on a large sample of photometrically selected 
galaxy clusters from the SDSS DR7. We detect a satel- 
lite alignment only when we use the isphotal PAs. As we 
noted in §3.3, the isophotal PA tends to trace the outer 
profile of the galaxy while the model fit PAs tend to trace 
the inner part of the galaxy. A direct interpretation of 
the measurement results could be that the outer part of 
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Figure 20. Fraction of projected field galaxies at different redshift 
bin of size 0.05. It maintains constant except in the lowest redshift 
bin. 



the satellite galaxy is more susceptible to the the gravi- 
tational torque and thus shows an orientation preference 
toward the BCG. However the inner part of the galaxy is 
not affected much by the tidal torque and does not show 
preference toward the BCG. The measured discrepancy 
of the satellite alignment from different PAs could be a 
manifestation of the twisting of galaxy shape from inner 
part to outer part. However, another possibility of this 
discrepancy could be that the light from BCG contami- 
nates the measurement of the PA based on the isophote 
fit to the outer region of the galaxy and lead to a "arti- 
ficial" alignment. By comparing the dependence of 5 on 
BCG apparent and absolute magnitudes, we favor the 
latter explanation. This means that, though the tidal 
torque within the galaxy cluster may induce the satellite 
alignment, we are not yet able to detect them based on 
our current SDSS data. It will be definitely an interest- 
ing question to address with the forthcoming high qua! 



ity data such as that from the Dark Ener gy Survey ( |The 
Dark Energy Survey Collaboration||2005 ). 



For the BCG alignment, by introducing the alignment 
parameter 7, we detect a strong redshift and BCG abso- 
lute magnitude dependences of the alignment. The red- 
shift dependence cannot be explained by our known sys- 
tematics. This result implies that the BCGs orientation 
is a dynamically evolving process and gets stronger as the 
cluster system evolves. For the dependence of 7 on the 
absolute magnitude of BCG, our result is qualitatively 
consistent with the conclusion that clusters with BCG 



dominance show st ronger BCG alignment in (Niederste- 
Ostholt et al.|2010 1 . Furthermore, based on a subsample 
of the BCGs whose stellar masses are available, we show 
that the BCG alignment signal becomes stronger as the 
BCG stellar mass increases. This result indicates that 
more massive BCGs (with lower absolute magnitude) are 
more likely to align with the major axes of clusters. 

We must take great caution when interpreting the de- 
pendence of 7 on BCG absolute magnitude and stellar 
mass since the purity of the cluster sample may also de- 
pend on the BCG absolute magnitude and stellar mass. 
As the cluster purity decreases, the alignment signal will 
decrease too. The faintest two bins in Figure|l4 show null 
alignment signal, which may also be due to the signifi- 
cantly decreased cluster purity. Nevertheless, we can still 
see a trend that 7 increases as the BCG absolute magni- 
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tude increases by looking at the bright end of the sample 
where we are confident about the cluster purity. Eval- 
uating the cluster purity variation w.r.t BCG absolute 
magnitude turns out to be difficult because it requires a 
mock galaxy catalog that has BCG information properly 
built in. The way the mock catalog is constructed will 
impact the results significantly. Therefore, we think the 
best way to check this purity variation w.r.t. magnitude 
is to perform similar analysis with deeper data in the 
near future, suc h as the data from the upcoming Dark 
Energy Survey (The Dark Energy Survey Collaboration] 



2005). 



REFERENCES 

Abazajian, K., & Sloan Digital Sky Survey, f. t. 2008, ArXiv 
e-prints 

Agustsson, I., & Brainerd, T. G. 2006, ApJ, 644, L25 
Austin, T. B., & Peach, J. V. 1974, MNRAS, 168, 591 
Azzaro, M., Patiri, S. G., Prada, F., & Zentner, A. R. 2007, 

MNRAS, 376, L43 
Bernstein, G. M., & Jarvis, M. 2002, AJ, 123, 583 
Bernstein, G. M., & Norberg, P. 2002, AJ, 124, 733 
Binggeli, B. 1982, A&A, 107, 338 
Brainerd, T. G. 2005, ApJ, 628, L101 
Brown, F. G. 1938, MNRAS, 98, 218 
Carter, D., & Metcalfe, N. 1980, MNRAS, 191, 325 
Catelan, P., Kamionkowski, M., & Blandford, R. D. 2001, 

MNRAS, 320, L7 
Ciotti, L., & Dutta, S. N. 1994, MNRAS, 270, 390 
Ciotti, L., & Giampieri, G. 1998, ArXiv Astrophysics e-prints 
Colless, M., et al. 2001, MNRAS, 328, 1039 
Croom, S. M., et al. 2009, MNRAS, 392, 19 
Dressier, A. 1978, ApJ, 226, 55 

Faltenbacher, A., Jing, Y. P., Li, C, Mao, S., Mo, H. J., Pasquali, 
A., & van den Bosch, F. C. 2008, ApJ, 675, 146 

Faltenbacher, A., Li, C, Mao, S., van den Bosch, F. C, Yang, X., 
Jing, Y. P., Pasquali, A., & Mo, H. J. 2007, ApJ, 662, L71 

Faltenbacher, A., Li, C, White, S. D. M., Jing, Y., Shu-DeMao, & 
Wang, J. 2009, Research in Astronomy and Astrophysics, 9, 41 

Hao, J., et al. 2009, ApJ, 702, 745 

— . 2010, ArXiv e-prints 

Hawley, D. L., & Peebles, P. J. E. 1975, AJ, 80, 477 

Hirata, C. M., Mandelbaum, R., Ishak, M., Seljak, U., Nichol, R., 

Pimbblet, K. A., Ross, N. P., & Wake, D. 2007, MNRAS, 381, 

1197 



Kauffmann, G., et al. 2003, MNRAS, 341, 33 
Kim, R. S. J., et al. 2002, AJ, 123, 20 

Knebe, A., Gill, S. P. D., Gibson, B. K., Lewis, G. F., Ibata, 

R. A., & Dopita, M. A. 2004, ApJ, 603, 7 
Kuhlen, M., Diemand, J., & Madau, P. 2007, ApJ, 671, 1135 
Mandelbaum, R., Hirata, C. M., Ishak, M., Seljak, U., & 

Brinkmann, J. 2006, MNRAS, 367, 611 
Niederste-Ostholt, M., Strauss, M. A., Dong, F., Koester, B. P., 

k, McKay, T. A. 2010, MNRAS, 405, 2023 
Paz, D. J., Sgro, M. A., Merchan, M., & Padilla, N. 2011, 

MNRAS, 502 

Paz, D. J., Stasyszyn, F., & Padilla, N. D. 2008, MNRAS, 389, 
1127 

Pereira, M. J., & Bryan, G. L. 2010, ApJ, 721, 939 

Pereira, M. J., Bryan, G. L., & Gill, S. P. D. 2008, ApJ, 672, 825 

Pereira, M. J., & Kuhn, J. R. 2005, ApJ, 627, L21 

Plionis, M., Benoist, C, Maurogordato, S., Ferrari, C, & 

Basilakos, S. 2003, ApJ, 594, 144 
Salim, S., et al. 2007, ApJS, 173, 267 
Sastry, G. N. 1968, PASP, 80, 252 

Siverd, R. J., Ryden, B. S., & Gaudi, B. S. 2009, ArXiv e-prints 

Stoughton, C, et al. 2002, AJ, 123, 485 

Strauss, M. A., et al. 2002, AJ, 124, 1810 

Struble, M. F., & Peebles, P. J. E. 1985, AJ, 90, 582 

The Dark Energy Survey Collaboration. 2005, ArXiv 

Astrophysics e-prints 
Tormen, G. 1997, MNRAS, 290, 411 
Tucker, G. S., & Peterson, J. B. 1988, AJ, 95, 298 
Ulmer, M. P., McMillan, S. L. W., & Kowalski, M. P. 1989, ApJ, 

338, 711 

Vitvitska, M., Klypin, A. A., Kravtsov, A. V., Wechsler, R. H., 
Primack, J. R., & Bullock, J. S. 2002, ApJ, 581, 799 

Wang, H. Y., Jing, Y. P., Mao, S., & Kang, X. 2005, MNRAS, 
364, 424 

Wang, Y., Park, C, Yang, X., Choi, Y., & Chen, X. 2009, ApJ, 
703, 951 

Wang, Y., Yang, X., Mo, H. J., Li, C, van den Bosch, F. C, Fan, 

Z., & Chen, X. 2008, MNRAS, 385, 1511 
Yang, X., van den Bosch, F. C, Mo, H. J., Mao, S., Kang, X., 

Weinmann, S. M., Guo, Y., & Jing, Y. P. 2006, MNRAS, 369, 

1293 

York, D. G., et al. 2000, AJ, 120, 1579 

Zentner, A. R., Kravtsov, A. V., Gnedin, O. Y., k, Klypin, A. A. 
2005, ApJ, 629, 219 



ACKNOWLEDGMENTS 

JH thanks Scott Dodleson for helpful comments and Eric Switzer for helpful conversations. Funding for the SDSS 
and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Sci- 
ence Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese 
Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web 
Site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participat- 
ing Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, 
University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the 
Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear 
Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese 
Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), 
the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University 
of Washington. 



